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Abstract 
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r~| ' In the limit of a nonlinear diffusion model involving the fractional Laplacian we get a 

■ "mean field" equation arising in superconductivity and superfluidity. For this equation, we 
obtain uniqueness, universal bounds and regularity results. We also show that solutions 
with finite second moment and radial solutions admit an asymptotic large time limiting 

■ profile which is a special self-similar solution: the "elementary vortex patch" . 
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1 Introduction 

a • We consider the class of evolution models with nonlocal effects given by the system 

(1.1) Ut = V ■ {uVp), p = ICu. 

where u is a function of the variables (x, t) to be thought of as a density or concentration, and 
therefore nonnegative, while p is the "pressure" and v = —Vp represent the mass propagation 
speed. We postulate that p is related to u via a linear integral operator /C which we assume in 
practice to be the inverse of a fractional Laplacian, that is, JC = (— A)"'^, with s G (0, 1]. More 
general integral operators are also now being considered in the literature. 

The problem is posed for t > and x G M"", n >2. (The case n = 1 is easier and a bit different, 
and a short explanation will be given below). Finally, we give initial conditions 

(1.2) u{x,0) = uo{x), a; G M", 
where Uq is a nonnegative integrable function in M". 
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Nonlocal diffusive models. These correspond to the equation 

(1.3) = V- (nV(-A)-%), 

in the parameter range < s < 1. This equation can be considered as a nonlocal version of the 
standard porous medium equation ut = A(n^). For convenience, we will often refer to it as the 
FPME (fractional porous medium equation). 

The FPME has been derived starting from the continuity equation 

(1.4) dtu + V ■ (uv) = 0, 

where m is a density and v is a velocity field, which according to Darcy's law derives from a 
pressure, v = — Vp (as in the theory of gases in porous media). Physical applications range from 
macroscopic evolution of particle systems with short- and long-range interactions [3B], to phase 
segregation [31], [32], and dislocation dynamics [35], [15]. The existence of a weak solution, the 
relevant integral estimates, as well as the property of compact support have been established by 
Caffarelli and the second author in the paper [20] . The existence is shown by passing to the limit 
in the sequence of smooth solutions of approximate problems, and using the compactness of the 
integral operators involved in the weak formulation and its energy estimates for s < 1 (thus the 
method does not directly apply to s = 1). Note that no uniqueness is proven, and examples are 
constructed where the maximum principle does not hold for s > 1/2. 

Hydrodynamic or mean field equation. The case s = 1 of equation f ll.ll) - ([L2l ) 

(1.5) Ut = V-{uVp), p={-A)-^u, 
supplemented with initial data: 

(1.6) u{x,0) =u°{x) , 

is no more a diffusion equation. This equation was first studied in by Lin-Zhang [H]. There, 
existence is proven by a vortex point approximation, and uniqueness by estimates in Zygmund 
spaces. In dimension 2, the equation is directly related to the Chapman- Rubinstein-Schatzman 
mean field model of superconductivity [23] and to E's model of superfluidity [26], which would 
correspond rather to the equation 

(1.7) ut = V-{\u\Vp) 

which coincides with (11.51) when u > 0. That model is a mean field model for the motion 
of vortices in a superconductor in the Ginzburg-Landau theory. There, u represents the local 
vortex-density, and p represents the induced magnetic field in the sample. In this application 
the equation is really posed in a bounded domain and for signed "vorticity measures" u (which 
can be measures instead of functions), but the study of fll.5p in with u > 0, corresponding to 
the situation where the vorticity is positive, is a reasonable start. It can also be viewed as the 
gradient-flow version of the Euler equation in vorticity form i.e. 



ut = V-{uV^p) p={-Ay\, 
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hence also the term hydro dynamic equation. In |2], Ambrosio and the first author studied 
the model in a bounded two-dimensional domain (taking into account the possibility of flux of 
vortices through the boundary), via a gradient flow approach: the equation is the gradient-flow 
for a natural energy (see below) for the 2-Wasserstein distance. There, u is denoted /i, and p is 
denoted h. Solutions with u measure (and not L^) are also treated there. The equation (11.71) with 
changing sign u has been studied by Masmoudi- Zhang ^9] in M^, Ambrosio-Mainini-Serfaty [3] 
in bounded domains of the plane and in the whole plane, again using a gradient-flow approach. 

This equation is also related to nonlocal aggregation models, which have attracted a lot of 
attention lately, notably in conjunction with the Keller-Segel model [37], see e.g. (TUl [HI [12] • 
They correspond to (11. ip with K, a general integral operator corresponding to an attractive 
(instead of repulsive) interaction. Particularly relevant to us is the paper of Bertozzi-Laurent- 
Leger [12] which deals mostly with aggregation via the Coulomb kernel, but contains a section 
(Section 3) that deals with the repulsive case, i.e. exactly the equation (ll.5p -f [L6|) . In addition 
to existence and uniqueness results, it makes a detailed analysis of the property of support 
propagation and asymptotic behaviour for solutions with compactly supported initial data. 

Outline of results. In the present paper we are interested in obtaining Equation ( 11. Sp as the 
limit of the nonlocal diffusive equations (ll.3p when s — )■ 1; and in expanding the theory for the 
limiting equation. We consider the initial- value problem posed in M", for all dimensions n > 2, 
with comments on the easier case n = 1. We develop an existence theory for weak solutions with 
a nonnegative finite measure as initial data. Such solutions are uniformly bounded for positive 
times in the form \u{x,t)\ <l/t. We prove a uniqueness result for bounded solutions. 

As another important feature of the theory, we derive the support propagation properties (for 
general solutions), and compare the result with a similar property for the fractional diffusion 
approximation . 

Some insight is obtained with the introduction of the special formulation of the equation under 
conditions of radial symmetry, leading to a non-homogeneous Burgers' equation for the so-called 
mass function. With this Burgers' equation approach, we recover the existence of a fundamental 
solution in the form of a round vortex patch of constant (in space) density, and show that they 
are attractors for the evolution of all radial solutions. We also find other explicit radial solutions 
that allow to derive a new non-comparison result. 

We also study the question of large time behaviour for another class of initial data where the 
restriction of compact support is dropped, more precisely integrable data with finite second 
moment. The method, which is interesting for its own sake, is an energy method that follows 
the now well-known entropy - entropy dissipation approach, see [H] and references. 

Plan of the paper. The paper is organized as follows: Section [2] is devoted to the definition 
of solutions for s < 1 and s = 1 and to recalling some of their properties (conserved quantities, 
etc). In Section [3 we show the existence of solutions to (11.50 by taking the limit s — )■ 1 in 
solutions to (II. ip . Section [His devoted to the proof of uniqueness for bounded solutions. Section 
[5] contains the universal decay estimate for the solutions. In Section [6] we investigate the theory 
in the class of radial solutions via the Burgers' equation approach. In Section [7| we show that the 
compact support property is preserved during the evolution, see also [12]. Section [S] discusses the 
question of large time behaviour for two classes of initial data where the restriction of compact 
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support is dropped: (i) integrable data with finite second moment, (ii) all solutions that are 
radially symmetric in the space variable. 

Open problems. Here is list of the main problems that remain open: obtaining uniqueness 
results beyond the class of bounded solutions, analyzing the long-term behaviour of all solutions 
with integrable data, and finally treating the case of sign changing solutions. 

Acknowledgments. We would like to thank J. A. Carrillo for pointing out to us the reference 
[12] after a first draft of this paper was written. S. S. thanks the Universidad Autonoma de 
Madrid for its hospitality that allowed this work to be completed. S. S. was supported by a 
EURYI award, and JLV by Spanish Project MTM2008-06326-C02-01. 



2 Definitions and first properties of the solutions 

It will be convenient to write for any s G (0, 1] the equation in the form 

(2.1) dtU = V ■ iuWK.u), ]Cs = {-A)-\ 

K,s is a positive essentially self-adjoint operator. For convenience, we introduce "Hs = K}J'^ = 
(— A)~'^/^. We write p = ICgU, so that u = {—AYp and HgU = {—Ay^'^p. 

We start with the definition of weak solutions. The notation A4~^{W^) denotes the space of 
positive Radon measures on M" with finite mass. 

Definition 2.1. For any s G (0, 1], a weak solution of (12. ip in Qt = K" x (0, T) is a nonnegative 
function u{x,t) such that u G L^{Qt), p = ICsU E L\{t,T),WI;^{W)), uWp E L\{t,T) x M") 
for all T > 0, and the identity 

(2.2) j j {u(pt-u\/p-W(p)dxdt = Q 

holds for all continuously differentiable test functions compactly supported in Q = M" x (0,T). 
If s < 1 we also assume that u is continuous in Qt- 

Definition 2.2. For any s G (0, 1], a weak solution of the problem formed by equation (12. ip in 
Qt = M" X (0,T) with initial data fi G A^^(M") is a continuous nonnegative function u{x,t) 
with u G L\Qt), p = lCsUe L^iQ.T),^^!^^^)), u\/p e L^{Qt), and the identity 

(2.3) j j {u(j)t-uVp- dxdt + J 0(x, 0) dfi{x) = 

holds for all continuously differentiable test functions m M" x [0, oo) such that cj) has compact 
support in the space variable and vanishes for t > T. 

In practice the solutions obtained below have much better properties that can be exploited in 
the analysis. In particular they will be uniformly bounded for t > r > and satisfy energy 
inequalities. Following |33] we may call such improved solutions dissipative solutions. 
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2.1 Properties of the solutions in the case s < 1 

The existence of weak solutions in the sense of Definition 12.21 is proved in [20] under the as- 
sumption that the initial data is actually a bounded function Uq > decaying exponentially 
as |x| — )■ oo, cf. Theorem 4.1 of [20]. It is extended to general nonnegative initial data in 
[T9] . Let us list the main properties of those solutions with attention to the way they depend 
on s G (0, 1). We point out that we do not use all of them in the sequel. 

- Conservation of mass: The solutions exist globally in time and for every t > we have 
(2.4) J u{x,t)dx = J UQ{x)dx. 

- Conservation of positivity: uq > implies that u{-,t) > for all times. 

-L^ ESTIMATES: the norm of the constructed solutions does not increase in time, for any 
1 < p < oo. 

- A general comparison theorem, i.e., a form of the usual maximum principle, is proven to 
be false by constructing counterexamples for s > 1/2, see [20] • However, comparison of solutions 
for an integrated version of the equation is established in dimension = 1 by the method of 
viscosity solutions in Biler-Karch-Monneau [15]. We will return to that issue below. 

- Finite speed of propagation. If the initial function has compact support so does the 
solution at any given time t > 0. Estimates on the growth of the support for very large time are 
given in [20] and [21] . 

-Log-entropy estimate: We have 

therefore, solutions with uologuo G L^(M") stay in the same space and moreover satisfy 

I VHomP dxdt < oo. 



- Energy estimate. We first make the observation that if p = KgU we have 

(2.6) E(u)-.^ju,,.^j((-^)"yf,.^j(n.uf,.>,. 

This formal computation is justified for the class of constructed weak solutions. This is the main 
energy we are going to use. We also have formally 

/|VH.«P..^/|V(-Ar/W..^/.(-Ap)... 

The estimate on the evolution of the energy according to our flow is then 

— I \'Hsu\'^ dx = 2 I {T-Lsu){Hsu)tdx = 2 j {lCsu)utdx = —2 f u\\/}Csu\^dx. 
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where the last identity is obtained after using the equation and integrating by parts. This is 
justified for smooth solutions, while in the generality of weak solutions we get 



(2.7) [ \'Hsu{t2)\^ dx + 2 I I u\VlCsu\^ dx < f \'Hsu{ti)\^ dx 

- BOUNDEDNESS: The solutions to f l2.ip constructed in [20] have the following boundedness 
effect (used in the subsequent paper [IH]): For all < s < 1, weak energy solutions in are 
bounded for all positive times. More precisely, there is a constant C = C{n,s) such that 

(2.8) ||n(t)||oo<C(n,s)t-"||no||Ii 

Dimensional considerations imply that a = n/{n + 2 — 2s) and cr = (2 — 2s) /{n + 2 — 2s). This 
fact is proved in [19]. Notice that as s — 1 we have a — )■ 1 and cr — )■ 0. Once we can control the 
constant C{n, s) uniformly in s we can pass to the limit s — > 1 to obtain the a priori estimate 
(already found by Lin and Zhang in their study in dimension n = 2 in ^^) 

(2.9) u{x,t)<C/t 

for the solutions of the limit equation obtained in the limit process, this will be done in the next 
section in Theorem 13.11 This just requires a closer inspection of the proof of [T9] . 

- Energy solutions: The basis of the boundedness analysis is a property that goes beyond 
the definition of weak solution. We will review the formulas with attention to the constants 
that appear since this is not done in [19]. The general energy property is as follows: for any F 
smooth and such that f = F' is bounded and nonnegative, we have for every < ti < t2 < T, 

j F{u{t2))dx- J F{u{ti))dx =-J ' J V[f{u)]uVpdxdt = 

where /i is a function satisfying h'{u) = u f'{u). We can write the last integral as a bilinear form 

(2.10) j Vh{u)V{-A)-'u dx = Bs{h{u),u) 
and this bilinear form Bs is defined on the Sobolev space iy^'^(M") by 

(2.11) B,{v,w) = Cn,s jj ^v{x) _'^^^^,^_^ Vw{y) dxdy = j j M-s{x,y)Vv{x)Vw{y) dx dy 

where J\f-s{x,y) = Cn,s\x — y\~^"'~'^^^ is the kernel of operator (— A)~*. After some integrations 
by parts we also have 

(2.12) Bs{v, w) = Cn,i^s J I {v{x) - v{y)) ^^_J^,,^2(i-s) (^(^) - ^(2/)) 

since — AA/l^ = -^i-s- It is known [SH] that Bs{u,u) is an equivalent norm for the fractional 
Sobolev space W^~'^''^(W') and we will see just below that CnA^s ~ Kn{^ — s) as s — > 1, for some 
constant Kn depending only on n. 
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Proposition 2.1. The constant C{n,s) in fl2.8p is uniformly bounded as s —> 1, i.e. for any 
solution to fl2.ip . we have 



Ht)\\oo < 



Cnll'^olljl 



fa 

with a = n/{n + 2 -2s) and a = (2 - 2s) /{n + 2 -2s). 

Proof. It is based on a careful scrutiny of the ideas and proof of paper [19], Section 4. 

(i) We have Cn,i-s ~ Kn (1 — s) as s — )■ 1. Indeed, if < s < 1 we can also use the integral 

representation 



(2.13) {-^r9{x) = Cn,s P.V. / f''^ dz = K* g( 



x] 



where P.V. stands for principal value and Cn,s = ^J^J/i^p^'fly^ is a normalization constant, see for 
example [SI [56]. Note that Cn,s ~ KnS as s — and Cn,s ~ C'^{1 — s) as s — 1. 



(ii) Next, we need the following result from [50 



lim r / / (^(^) ""^y^^' dxdy = I u{x)' dx. 



r^0+ J J \x — 

Setting r = 1 — s and combining this with fl2.12p and Cn,i-s ~ Kn{^ — s) as s — 1~, we deduce 
that 

lim Bsiu,u) = K„|§"^^|||m||2. 

s->l- 

(iii) We can now revisit the proof of boundedness of Theorem 4.1, in Section 4 of [19], and check 
that the constants in all the arguments are uniform in s as s — )■ 1. □ 

We end this section with some further properties of the solutions. 

- Holder estimates are proven for bounded energy solutions in [19] and hold for all s G (0, 1). 
The Holder exponent has to deteriorate as s — 1 since we do not expect solutions of the limit 
equation to be continuous, as can be seen in the case of the "vortex-patch solutions" below. 

- More general data. The property of boundedness and the regularity are used in [19] 
to obtain solutions with similar properties for all nonnegative initial data in L^(M"'). We will 
generalize that result to measures in Section |3] below. 



3 Existence of solutions as limits of fractional diffusion 

In dimension n = 2, existence and uniqueness of positive L°° solutions has been proven by Lin 
and Zhang in [H] (they also proved existence of positive measure- valued solutions). Existence 
with positive initial data of finite energy has also been proven by a gradient flow approach, in 
bounded domains of the plane in [2], and in all of in [3]. In general dimension, existence can 
be obtained by taking the limit s — )■ 1 in the solutions for s < 1 constructed in [20] as we will 
see below. 
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We examine here the situation s = 1 as hmit of the equations for s < 1. Since all the estimates 
above are uniform in s as s — )■ 1, we may write the formal version for s = 1, where the equation 
is fll.Sp . Here we need to distinguish the case n > 3 from the case n = 2 which we will examine 
separately below. For n > 3, and noting that now V,u = {-Ay/^p, the energy is 

(3.1) E{u)= / updx= {nu)'^dx= / \Vp\'^dx, 

which is an analogue of the energy of [2] (there, it was set in bounded domains) for which the 
equation is a Wasserstein gradient flow, and its evolution is given by 

(3.2) ^ = -2 [ u\Vp\^dx. 

dt J^n 

This is a formal computation that will be justified as a consequence of the passage to the limit 
s — 1 in formula ( 12. 7p . This energy would not make sense for n = 2 because, unless J u = 0, 
p = {—A)~^u is never of gradient in in dimension 2 (the Green kernel is a logarithm), however 
a definition is given below in Section 13.21 



3.1 The existence result by s — > 1 limit 

We next prove the following general result in all dimensions n > 2. 

Theorem 3.1. For every fi G Ai^(W^), there exists a continuous and nonnegative weak solution 
of equation (11.51) . with data /i in the sense of Definition \2.S\ It can be obtained as the limit s — > 1 
of the solutions of the FPME with < s < 1 with data e L^(]R") fl L°°{W) that approximate 
/i in the sense of measures. In addition these solutions satisfy 

(3.3) |k(t)|U°c<^ 
for some C > depending only on n. 

Proof. Existence is proved by passing to the limit s — )■ 1 in the solutions of the FPME (11. ip . 
The outline is similar to the proof of Theorem 13.1 in [IS] where s was considered constant, but 
we cannot use the continuity in time of the solution, since it may be lost in the limit s — > 1. 

Let us consider initial data m° = /i in A^+(]R"), and a sequence of nonnegative functions 
approaching as measures as /c — )■ oo. (If m° G we may simply take = m°.) Let us then 
consider the family of solutions Uk{x, t) with parameter — 1 of (11. ip given by Theorem 4.1 of 
[20] or Theorem 13.1 of [19], and let pk{x,t) = {—A)~'''uk{x,t) be the corresponding pressures. 

(i) In a first step, we assume that n > 3, and prove that the limit exists (along subsequences 
Skj — > 1) and is a solution of the equation in the sense of Definition 12. H without discussing the 
question of the initial data. By Proposition 12.11 we have that for any r > 0, 

(3-4) lkfe||L°°((0,T),Li) < C 
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where C > only depends on r. In order to pass to the hmit in the weak formulation 

(3.5) Jj {uk (l)t - UkVpk ■ Vcj)) dxdt = 

for compactly supported in x (0,T), it suffices to have weak convergence of Uk, together 
with strong convergence of Pk and Vpk- Let us work on [r, T) containing the support of 0. From 
(13 ■4p we have weak convergence (up to subsequence extraction) of Uk in L'^([r, T) x Bfj) for any 
q>l, R>0. 

- For Pk, we may invoke the embedding theorems associated to the Riesz kernels, which tell us 
that Pk = {—A)~^Uk is uniformly bounded in L°°((r, T), LI^Br)) for all large q and i? > 0. More 
precisely, q> n/{n — 2s), which is uniform in s ~ 1 if n > 3. 

- Moreover, Wpk has a similar bound in spaces of the form L°°((t,T), W[^f) with s' < 2s — 1 
and convenient q > n/{n — s). Since we are interested in s — 1, we can restrict to 2s — 1 > 
so we can take s' > and both [pk] and (Vpfc) are a locally compact family with respect to the 
space variable. 

- We still need compactness in time for pk and Vpfc, and this is obtained from Aubin-Simon's 
results, cf. [B], [SS]- First, for p^, we use the facts that 

(3.6) Pk is bounded in L°°((r, T), U'{Br)) 

dtPk = {-Ay^'^dtUk = (-A)-^'= V ■ (ukPk) 
and since UkPk is bounded in L'^{Br) for any g > and uniformly in t G [r, +oo), it follows that 

(3.7) dtPk is bounded in L°°((r, T), W'^'^-^'\Br)). 

(13. 6p and (13.71) ensure, via Aubin-Simon's theorem, that pk is compact in C{\t,T], L'^{Bpi)), for 
any g > 1 and i? > 0. The proof for Vpfc is entirely similar: first, from the above, we have 

(3.8) Vpk is bounded in L°°((r, T), W''^'\Bn)) 

dtVpk = {-AydtUk = V(-A)-^^V ■ {ukPk) 
and since UkPk is bounded in L'^{Br) for any g > and uniformly in t G [r, T], it follows that 

dtVpk is bounded in L~((r, T), W'^"'-'^'\Br)). 

Again, these suffice to apply Aubin-Simon's result and obtain that Vpk is compact in say 
C{\t,T], L'^{Bji)) for every r > 0, T > r, i? > 0, g > 1. The limits after extraction of Uk 
and Pk are denoted u and p. 

Taking smooth compactly supported in space and vanishing for t > T, it follows from these 
convergence results that 

(3.9) hm / / {uk4>t — UkVpk ■ V0) dxdt = / {u(j)t — uVp ■ V0) dxdt. 
I'^'^Jo J Jo J 
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(ii) When n = 2 we do not have a uniform bound of the pressures pk in L'^ spaces vahd for all 
s ~ 1, since this is a limit case of the Sobolev embeddings. Fortunately, the pressures do not 
appear directly in the equation, but only through the pressure gradients for which the estimates 
are uniform in the range s f^:^ 1. The rest of the argument follows. 

(iii) There remains to control terms corresponding to the time interval [0, r] in order to complete 
the proof that the solution takes on the initial data, as in Definition 12.21 For that we follow |19] . 
proof of Theorem 13.1. We consider (p G C^(W^ x [0,oo)) compactly supported in space and 
vanishing for t >T. First, let r > be arbitrary. Applying the result of (13. 9p we have that 

J {uk4>t — Uk'Vpk ■ V0) dxdt = J J {u(f)t — uVp ■ V0) dxdt. 

In addition, by weak convergence of to we have 

(3.11) hm j ul{x)(f){x,0)dx = J J u^{x)4>{x,0) dx. 

There remains to show that the contributions on the interval [0, r] can be made small. 
Starting again from Proposition 12.11 we have 

\\uk{t)\\Loo < C(Mo)r" a = n/{n + 2 - 2sk), 

which yields by conservation of mass 

WMmi. < \\uk{t)\\L4ukmi-J < c{uo)r^^p^'\ 

By Riesz embedding we deduce 

/ \vnsMdx < c\\u,{t)\\l, < cr°(^-^^i, ^ = - - - 

J 2 p n 

where we recall HsU = {—A)~'^^'^u for any m, and combining the definitions of a and p, this yields 



Combining with the energy inequality (12. 7p . we deduce that we also have 
(3.12) [ [ u,\Vpk\^ < Ctm^-\ 



We may now bound the remaining contribution JJ" J {uk(j)t — u^Vpk ■ dxdt as follows: for 
tm = 2"*", we first write 



^m — 1 



Uk(pt 
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and second 



2 / rtm-l f \ 2 



1 



^ j."+2-2si. 2 ^ 9/^+27. 2 n+2-2sj, n+2-2sj. 

— m— 1 — ^mi^m — ^ i^m 



Summing over m we find that the contributions on [0, r] tend to as r tends to and combining 
with dXTUD - flXTT]) we obtain 



00 







hm / I {uk4>t — pk ■ <p) dxdt + / Uy.{x)(t){x,Q) dx 



as desired. 



(iii) Finally, by passing to the limit — > 1 in fl2.8p and using the fact that C {n, s) does not 
blow up in this hmit, we obtain the universal estimate ( I3.3p . □ 



Remark 3.1. The handling of the initial data in the proof above could be replaced by showing 
that the initial data is taken via a result of continuity of the solution in Wasserstein metric, see 
Lemma \5. 2\ below. 

When Sfc is taken fixed this is the proof of existence of solutions with measure data that improves 
Theorem 10.1 of [19]. 



Theorem 3.2. For every fi G A^^(M"), there exists a continuous and nonnegative weak solution 
of FPME, < s < 1, with data /i in the sense of Definition\KE It is a bounded energy solution 
for allt>T>0. 

Remark 3.2. We will provide in Proposition I5.il below another approach to the existence of 
solutions for s = 1 via the vanishing viscosity method. 



3.2 Definition of the energy in two space dimensions 

In dimension 2, even though p = [—A)~^u is not in general of finite energy J |Vpp < 00, we 
may still define an energy by a "renormalization" procedure. Here we reproduce the definition 
of [3], which itself follows a similar construction for Ginzburg-Landau in [51] and [T3] . 

The first idea is to define the energy in a "renormalized way" , over functions u such that j u = k 
with k a fixed constant (we may do this without loss of generality since the integral is preserved 
by the evolution) by 

(3.13) S{u) = hm ( [ \Vp\' -^hgn) . 

\JBiO,R) / 

Contrarily to / |Vpp, it is possible to make this energy finite even if /c 7^ 0. 
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Another definition is the following: consider any smooth positive compactly supported Uq with 
/]R2 ""o = k, and let po = {—A)~^Uo. For any u with J u = k such that w := p — po is in if^(M^), 
we define the energy by 

(3.14) E{u)= I \Vw\^ + 2 I uqW 

In fact, expanding w as p — po and noting that J^^^q |VpoP ~ l^^°§-^ + '^^^ as i? — oo, and 
using an integration by parts we see that if w G if^(M^), the two definitions differ only by the 
constant (depending only on the choice of Uq). So, for energy we may take E. 

Formally, along the flow, since Wt = Pt, "we have 

^E{u{t)) = 2 [ Vw-Vwt + 2 [ uowt 



wApt + 2 uopt 

2 Jr2 

w + po)ut = 2 put 
\Vp\^u. 

An alternate approach would have been to view the (— A)~^ operator as a limit of (—A + hl)~^ 
as /i — )■ 0. 



4 Uniqueness of bounded solutions for s = 1 

In this section we focus on obtaining a uniqueness result for s = 1 for bounded solutions. 
Uniqueness is not expected in all situations, see [1]. We also recall that there is no maximum 
principle so uniqueness cannot follow from comparison arguments. The proof works only for 
s = 1; no uniqueness proof is known for the FPME in dimensions n >2. 

A natural space in which to work for such continuity equations is the space of probability mea- 
sures endowed with the Kantorovich- Rubinstein- Wasserstein distance (commonly called Wasser- 
stein), indeed the PDE is a gradient flow for the 2- Wasserstein distance (point of view which 
was exploited in [21 |3]. For completeness, we recall (see e.g. [1]) that the p- Wasserstein dis- 
tance between two probability measures fi and u on M" (and by extension between two positive 
measures with same mass) is defined by 

(4.1) Wp{fi,iyr= mi J[ \x - y\P dj{x,y) 

where T{fi, v) is the set of probability measures on x M" whose first marginal is /i and whose 
second marginal is v. This inf is achieved by an "optimal transport plan" , and when 1 < p < oo, 
Wp defines a distance. 
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Theorem 4.1. There exists at most a unique solution of the initial-value problem for Equation 
ffT3]) m L°°((0,T) : L°°(M")). 

When n = 2, this theorem recovers with a different method the resuh on existence and unique- 
ness of bounded solutions of [H] but without requiring initial compact support. 



Proof. The proof is inspired by the method of Loeper [17] (see also [2]). We will need the 
following result, adapted from ^8] Lemma 8.1]. The result there is only stated in dimension 
n = 2, we include here a proof in any dimension n. 

Lemma 4.1. Let p = {—A)~^u. There exists a constant C > depending only on the dimension 
n such that 

Vx, y G M", \Vp{x) - Vp{y)\ < C(||no||Li + ||%||oo)|x - y\ (I log \x - y||l|,_,y|<i/e + l) 
Proof. Using the integral expression for A~^, we have, for any n > 1, 
(4.2) Vp{x) = Cn I K{x - z)u{z) dz K{x) = 

Let us denote d = \x — y\. First, we treat the case d < 1/e. We may write Vp(x) — Vp{y) = 
Cnih + I2 + h) with 



Ji := / {K{x-z)-K{y-z))u{z)dz, h ■= {K{x- z) - K{y - z))u{z) dz 

^\B{x,4) J B{x,4)\B{x,2d) 

{K{x — z) — K{y — z))u{z) dz. 

B{x,2d) 

For z G M"\-B(x,4), and since (i = |x — ?/| < 1/e we have \K{x — z) — K{y — z)\ < C\x — y\ for 
some constant C. Hence we may write < C|a; — 

For z G B{x,4:)\B{x, 2d), by a mean value argument, and since d = \x — y\, we have \K{x — 
z) — K{y — z)\ < C ■ We may thus write 

f dz 

Ihl ^ C\x — y\\\u\\L°° / -j = — I/IIIwIIlooI logd\ = C\\u\\l^\x — y\ \ log |x — 

J B{x,4)\B{x,2d) \Z — X]"^ 



For z G B{x,2d) we have \K{x - z) - K{y - z)\ < \K{x - z)\ + \K{y - z)\ = ^^z^ + 
and we conclude 

f dz dz 
\h\ < \\u\\l^ / 1 [— r + 1 r— r < C\\u\\L^d = C\x - y\\\u\\L^. 

J B{x,2d) 1-^ ^1 \y ^1 

Combining these estimates on Ji, I2 and we obtain the result in this case. 
If c? > 1/e, we split the integral into two parts: 

h -= f {K{x-z)-K{y-z))u{z)dz, h ■= [ {K{x - z) - K{y - z))u{z) dz. 

Jw^\B{x,2d) J B(x,2d) 

For the first part, we use again the fact that for z ^ B{x, 2d) we have \K{x — z) — K{y — z)\ < 
C < C\x — y\ since c? > 1/e. We thus find | Ji| < C\x — y\\\u\\Li. For z G B{x, 2d) we argue 
just as above for I3 and find I/2I < C| This concludes the proof. □ 
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We now give the proof of the theorem. Let us consider two bounded solutions Mi,M2 on [0, T] 
with initial data and let pi,P2 be the corresponding pressures. 

We consider the characteristics equation (for i = 1,2) 

dX 

(4.3) ^(t, x) = -Vpi{t, Xi{t, x)) Xi{0, x) = X. 

According to [T] (see also Corollary 3.3 in [17]), if Ui solve f ll.Sp . then we have 

(4.4) = (X,)#uO. 

We then define 



Q(t) = J u%x)\Xi{t,x) - X2{t,x)\'^ dx. 



We have that Q{0) = 0, and from (521), {ui{t) , U2{t)) < Q{t). 
Let us study the evolution of Q. Differentiating Q in time yields 

Q'{t) = -2 / u%x)iX,{t,x)-X2it,x)) ■ (Vpi(t,Xi(t,x)) - Vp2it,X2{t,x)))dx 



hence using the Cauchy-Schwarz inequality 

1/2 / /• \ 1/2 



(4.5) ^\Q'it)\ < (^j u^\X^-X2\^dx^ (^j v!'\Vpi{t,X,) - Vp2{t,X2)\^ dx 

= j Vpi(t, X,) - Vp2{t, X2) p dx 

On the other hand, (/ M°|Vpi(t,Xi) - Vp2(t, ^2)^^)^^^ < T^^'^ + T^^^ where 

T, = j u'\Vp,it,X,) - Vpiit,X2)\'dx 

and 

T2 = j u'^\Vp,it,X2) -Vp2{t,X2)\'dx. 



1/2 



Note first that from (|44 



(4.6) T2 = j U2{x,t)\Vpi{t,x) -Vp2{t,x)\^ dx. 

But Theorem 4.4. in [17] asserts that, if Ui and UJ2 are two positive measures of same finite total 
mass and with L°° densities, then = — (A)~^(jjj satisfy 

1 

(4.7) ||V^i - V^2|U2(Kn) < (maxdlwillioo, ||a;2||L°°))'' W2{uJi,U2) 
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Applying this to Ui(t) and U2(t), using the uniform L°° bound on these solutions, and inserting 
into (I4.6p . we find 

\T,\ < max(||ni(t)||^, \Mt)rjW^{u,{t),u,{t)) < CW^{u,{t),u,{t)). 
There remains to bound Ti. Using Lemma [4. 1^ and denoting M = \\u^\\ii, we may write 

|Ti| < C{M^ + ||ni(t)||L)( j n°(x)|Xi - X^^log^^ \X^ - X2\l\x,-x,\<i/edx 

Let us define f{x) = xlog^(a;) for < x < 1/e and f{x) = 1/e for x > 1/e. It is easy to check 
that / is concave, f{x) < xlog^ x, and moreover from the above (and the boundedness of ui) 

\Ti\<c{^j u'^{x)f{\Xi-X2\^)dx + j n°(x)|Xi -Xsl^dx) 

Using the concavity of /, by Jensen's inequality, we have 

j u\x)f{\X,-X2\')dx<Mf Q ^\X,-X2\'dx^ =Mf{^)<CMQ{t)\og'Q{t). 

We conclude that |Ti| < CQ{t){log'^ + 1). Combining with the bound for T2, recalling that 
1^2^(^1,1x2) < Q, and inserting all these bounds into fl4.5p . we find 



Q'{t)<\Q'{t)\<CQ{t)Jl + \og' ^ 



Q{ty 



where C depends on M and the uniform bounds on ui and U2. 

We may now conclude by a standard integration argument, and since J — , = log(logf + 

a/1 + log^ v) (using a change of variables u = logt^), that Q{0) = implies Q{t) = for t G [0, T], 
and since W2 {ui{t) , U2{t)) < Q{t) and W2 is a distance, the conclusion follows. □ 



5 The universal bound for s = 1 

We now show that the norm of solutions is decreasing in time. Indeed, the equation can be 
written as 

Ut = —u^ — Vu ■ Vp. 

which is the sum of absorption and transport with speed v = — Vp. This expression suggests a 
universal decay in time of the form 

(5.1) \\u{-Moo<\ 
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which we will establish below. This agrees with the limit s — > 1 of the smoothing effect of 
Proposition 12.11 for s < 1, cf. Theorem 13.11 What is surprising is that the bound is universal 
in the sense that it does not depend on any norm of the data. Note that in dimension 2, this 
bound appears in [H] . We have a slightly more precise result in terms of the L°° norm of the 
data. The subject was taken up in the recent paper [12], which is very close to the ideas we 
explain next in some detail for the reader's convenience. 

Theorem 5.1. For every weak solution u with hounded initial data uq we have the estimate 

1 



(5.2) u{x,t)< 



t + T 



where r = l/sup^Mo(x). Therefore these solutions coincide (via the uniqueness result for 
hounded solutions, Theorem \A.l\) with those of Theorem \3.1\ 

Remark 5.1. A similar result holds for the lower hound of the negative values when dealing with 
signed solutions. 

We start by studying a viscous approximation of the equation (which, by the way, is like the 
repulsive version of the Keller-Segel equation for chemotaxis [37t [36]). 

Proposition 5.1. Let G L°°(M") fl L-'^(M") with vP > 0. There exists a unique solution to 

n\ / ^tUe = eAUe + V ■ (UsVpe) Pe = {-A)~^Ue 

^^•^^ \n,(.,0) = uO 

Moreover is smooth, satisfies 

(5.4) <u,{x,t) <\\u^\\loo 

and converges locally in L'^(R'^ x [0, +oo) for any 1 < q < oo, to u, the unique solution to 

(USD -(HI. 

We note that this "vanishing viscosity approach" is an alternative way of obtaining existence 
of solutions to (fT3|) -( lL6|) . 

Proof. This natural "viscosity approximation" was previously considered in Schatzle-Styles [53], 
for a very similar equation, but in a bounded domain of M^. Because the proof is tailored for a 
bounded domain, it is not immediate to adapt it to our case. In order to prove the existence of 
solutions to (15. 3p . we set up a fixed point argument. 

Take T > and let 

if= it,(x,t) gL°°(M"x (0,T)), t;>0, and VtG [0,T], / v{x,t)dx< I u^{x)dx\. 
I Jr" Jr" J 

K is clearly a closed convex subset of the Banach space L°°(M" x [0,T]). 
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Lemma 5.1. Let v? E L°°(]R") nL^(M"). For any v E K, letting = (—A) ^v. Then, we have 
Vp^ e L°°(M" X [0,T]) and there exists a unique solution u := T{y) over [0,T] of the problem 



(5.5) 



dtu = eAu + V ■ {uVpv 
m(-,0) = 



r/iai solution satisfies T{v) G A' and ||T(f)||co.a <C.T is continuous and compact from K to 
itself. 

Proof. First we note that if v G i^, then for all t G [0,T], G fl which implies that 

Vp„ G L°°(]R"' X [0,T]) by Nash-Sobolev embedding, or more precisely, by observing 

\Vp^\{x,t)< I ^^^^^^rfy < ||t;||L^(Knx[o,T]) / „i + / v{y)dy 

JR" F ~ ?/| Jb{x,i) F ~ 2/1 ^M" 

< C (||t^||L°o(Rnx[0,r]) + ' 

where C depends only on the dimension. In fact, direct inspection allows to prove in addition 
that Vp^(-,t) G CO'"(M") for some a > 0. 

We are thus led to solving a linear equation of the form 



(5.6) 



dtu = eAu + V ■ (uf) 
m(-,0) = M° 



with / G L°°{MP X [0,T],M"). This can be solved by a standard convolution argument, see for 
example Lemma 12 in [16] which applies exactly (the fact that the dimension is n and not 2 does 
not affect at all the argument), and it provides us with a unique solution l~{v) > defined for 
all positive times, hence in [0, T]. The fact that for all t G [0, T], J^„ T{v) dx < /]g„ follows by 
integrating in time the equation. On order to show that T(f ) G K, there remains to show that 
T{v) G L°°(M" X [0,T]). 

Let us consider u = T{v) solving (15. 6p and let : M ^> R be a smooth function. Using (15.61) . 
integration by parts and Ap^ = —v, we have 

(5.7) f 0(u) = / <f)'{u) {eAu + V ■ (uVp,)) 
'it Jr" J 

= j <P"{u)\Vu\^ - j u(P"{u)Vp^ ■ Vu 
= j </>"(«) I Vwp - Jvp^- V[H{u)] 
= -e f d)"(u)\Vu\^ - I vHiu) , 



where we define H{s) = j^t(f)"{t) dt. If in addition is convex it follows that H > and, since 
f > 0, we find 

^ I m < 0. 
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Taking nonnegative convex and supported over [||m°||/^oo, +oo) we conclude that < u{x,t) < 
\\u^\\loo for every {x,t) G M" x [0,T]. This proves that T(f) G L°° and it belongs to K. T thus 
maps K into itself. It is also easy to see that it is continuous. 

Now, standard regularity theory (cf. |22]) gives us Schauder estimates 

||T(tO||cO.«(R"x[0,r]) < C {\\v\\L°°iR"x[0,T]) + IK' || Li(K" x [O.T])) < CT (||n°||ioo + . 

Since C°'° embeds compactly into L°°, T{K) is thus precompact. □ 



We may now apply Schauder's fixed point theorem |l3j which gives us the existence of a fixed 
point for Thence a solution to (15.31) . satisfying > and ||iie||L°°(iR"x[o,T]) < 

Finally once we have established (15. 4p . convergence of Ue to a solution of (II. 5p can be proven 
following the same steps as in the proof of Theorem 13.11 (except replacing Sk by 1). The solution 
is unique by Theorem 14. 1| which implies convergence of the whole sequence Wg. □ 

Proof of Theorem \5 . 1\ We first prove the desired estimate for solution to (15. 3p . which is 
smooth by standard parabolic regularity theory. Let us consider smooth such that 0(m) = 

for n < 0, (j)" > 0,0' > 0, and let us compute — / 0(ue — M[t))dx, with M some smooth 

nonnegative function of time. This yields 

(5.8) ^ / 0(m, - M{t))dx = [ ^'{ue - M{t))dtUe dx - M'{t) [ 0'(m, - M(t)) dx. 



dt „ 

Using the equation, we find, after several integrations by parts 



J (l)\Ue - M{t))dtUe = e J (P'iUe - M{t))AUe - J (p" {Ue - M{t))VUe " U^Vpe 

= -e j (t>'\u, - M(t))|Vn,|2 - j <f>"{u, - M{t)){u, - M{t))Vu, ■ Vp, 



- M{t) / 0"(n, - M{t))Vu, ■ Vpe 



< 



j V[H{u, - M{t))] ■ Vps - M{t) j <j)"{u, - M{t))Vu, ■ Vp,. 



Here, we have defined again H{s) = t(f)"(t) dt. Using —Ap^ = and another integration by 
parts, we find that the first term on the right-hand side is — / Hiu;; — M(t))us < 0. There 
remains to estimate in a similar way 

- M{t) J <P"{u, - M{t))Wue ■ Wpe = -M{t) j V[0'(n, - M{t))] ■ Vp, 

= -M{t) I <P'{u, - M{t))u, < -M{tf I <P'{u, - M{t)) , 



because 0' is supported on the positive part of the real line and M{t) > 0. We conclude that 
j 0'(n, - M{t))dtUedx < -M{tf j 0'(n, - M{t)) dx 
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and combining with fl5.8p 

(5.9) ^ j (j){ue - M{t))dx < -{M^{t) + M'{t)) j (t)'{ue - M{t)) dx. 

Choosing M{t) = ^ with r = (sup no(x))~^, we have M'{t) = —M'^{t), and it follows that 
— / (f){us - M(t)) dx < 0. On the other hand, since m° < M(0) we have / 0(m° - M(0)) dx < 0, 
and thus we have obtained that j (j){ue — M{t)) dx < for all t > 0. 

Taking now = 0„ to be a sequence that approximates the function (poix) = sgnj^{x). we 
obtain u^it) < M(t) a.e. for t > 0. Passing to the limit e — )■ 0, and using Proposition 15. H we 
arrive at the result. □ 

Remark 5.2. The estimate \u{x,t)\ < also applies to solutions of every sign which 

generalizes the case of "dissipative weak solutions" studied in l4^ . 

We conclude this section with a result of continuity of the solutions in Wasserstein distance, 
which we do not use however. 

Lemma 5.2. The solutions to ( \1.5\\ with vP G A^"^(]R") satisfying (13.31) are "absolutely contin- 
uous" (hence continuous) for the Wq-Wasserstein distance for any q G [l,oo).- more precisely 
Wq{u{t2),u{ti)) < jllm{s)ds with m G L\[0,T]) for any T > 0. 



Proof. By conservation of mass, for every t > we have the two a priori estimates ||M(t)||i < M, 
\W{t)\\oo < 1/t. By the Nash-Sobolev embeddings we then get Vp(-,t) G L°°(M"). A simple 
scaling argument then shows that 

(5.10) ||Vp(-,t)|U < C||Ap(-,t)|i;/"||Ap(-,t)i|(:-i)/« < CiMi/"t-("-i)/'^ 

The claimed continuity is a direct consequence of the nature of the equation (11.51) as a continuity 
equation, i.e. of the form dtU+V-{uv) = 0. Theorem 8.3.1 in [1] asserts that for such an equation, 
if \\v{t)\\Li{u{t)) dt < +00 then the solution u{t) is absolutely continuous in [0,T] with values 
in the space of measures equipped with the distance Wq. The precise definition of this (cf. [1]) 
is that Wq{u{ti),u{t2)) < m{s) ds with m G L^([0, T]), and it implies continuity in the usual 
sense. In the present case, v = and, from (I5.10p we have 

||v(t)|U.(,(t))cit = / (f |Vp(t)|%(t)) dxdt 
Jo \Jr" j 

< [ Ct^'"^-^ ( [ u{t) dx^ ' dt < C'Af^/^T^/". 
Jo \Jr" j 

It follows that the result of P applies and u(t) is continuous for the Wasserstein distance Wq. □ 
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6 Radial solutions. Burgers' equation 

In this section we analyze radially symmetric solutions and exhibit the special class of vortex- 
patch solutions. We use the Burgers' equation approach via the mass transform, which already 
appeared in [11] and has precedents in chemotaxis problems, [36]. The main purpose of most of 
this section is to fix the ideas and notation that we will use in the last two subsections, and in 
the asymptotic analysis of Section 18.21 

Lemma 6.1. Let u be a solution to (11.51) with radially symmetric bounded initial data uo(r), 
r = |a;|. Then u is radial, and the "mass function" defined by 

(6.1) M{r,t) = / s''-^u{s,t)ds 

Jo 

solves the Burgers' equation 

(6.2) Mt + ^MM, = , 

which after the change of variables s = r^/n can be written in the standard form Mt + M Mg = 0. 

Proof. This follows from elementary computations, or can be found in [TT]. □ 

Note that (16. 2p is now a local equation, contrary to ( 11. 5p and that it looks very much like the 
equation for chemotaxis, but with the opposite sign. The speed of movement is M/r"~^ = > 0. 
Notice that 

(6.3) Mt < 0. 

A common problem with Burgers' equations and other first order conservation laws is that they 
may develop shocks in finite time and the uniqueness of solutions needs the concept of entropy 
solutions, cf. [121 [38]. This is not the case here, due to the fact that M{r,t) is a monotone 
nondecreasing function of r and that the velocity v is positive, so that the characteristic lines 
spread. Conversely, it is easy to see that weak solutions M{s,t) of the Burgers equation with 
Ms > bounded, give rise to nonnegative, bounded, radial, weak solutions of the hydrodynamic 
equation ut = V{u'V{—A)^^u). 

Burgers' Equation (16.21) has the property of comparison (maximum principle) in the class of 
nonnegative and monotone solutions M. Since the equation is also 2-homogeneous in M, Benilan- 
Crandall's result [8] applies to give the pointwise inequahty 

(6.4) Mt > -y. 

Using the equation this implies the further inequality (that we already know) 

(6.5) Mr<—^, i.e. u{x,t)<-. 

These inequalities are optimal as we will show by means of the examples of the next section. 

The approach via the mass function M(r, t) will be used below to construct interesting classes 
of explicit solutions, that serve as prototypes or as counter-examples of the theory. 
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6.1 Solutions in dimension n = 1 

The Burgers' transformation, i.e., the formulation in terms of the mass function, apphes easily 
when we work in dimension n = 1 and we do not need the condition of radial symmetry to get 
in a direct way the standard Burgers' equation 

(6.6) Mt + MMr = 0, 

where we are still assuming that we deal with nonnegative solutions n > 0. 

6.2 The elementary vortex patch 

We next recover the existence of the vortex patch self-similar solution. 

Proposition 6.1. The equation (11. 5p admits the family of radial "vortex patch solutions" of the 
form 

(6.7) U{x,t) = j^X\x\<R{t+ry/" 

(as well as all their translates) where R > 0,t > are free parameters. 

We skip the proof of this known result, see for instance [12]. This configuration is called in the 
terminology of transport equations a rarefaction wave. The simplest form of the solution is 

(6.8) u{r,t) = jXr<m^/" 

which formally corresponds to r = 0. We call this solution the elementary vortex patch. Notice 
that it is a solution in the sense of Definition 12.21 It has a Dirac mass as initial data for the 
density, hence it qualifies as the fundamental solution of the problem. 

All of the above solutions can also be described as expanding mesas. The use of the term 
"mesa" appears in some limit cases to nonlinear diffusions, see e.g. [27] and [28] or [59] . 

Largest solution. Note that in the limit when the initial mass goes to infinity we obtain the 
solution 

(6.9) u{x,t) = J, M{r,t) = r'^nt 

that does not have an interface. This solution is an absolute a priori upper bound for all solutions 
of our problem with s = 1. It satisfies equations (16. 4p and (16.51) with equality. 

6.3 Self-similar solutions 

The vortex patch solution can be obtained as the limit of interesting self-similar solutions of the 
FPME with < s < 1. Such solutions are constructed in the paper [21] as self-similar solutions 
of the form 

(6.10) u{x,t) =t-''U{x/tl^) 
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In order to cancel the factors including t explicitly, we get the condition on the exponents 

(6.11) a + (2-2s)/3 = l 

If we also impose conservation of (finite) mass, which amounts to the condition a = n(3, we 
arrive at the precise value for the exponents: 

/3 = l/(n + 2-2s), a = n/{n + 2-2s). 

We also arrive at the nonlinear, nonlocal elliptic equation for the self-similar profiles 

(6.12) Vy ■ (t/Vy(P + a|?/n) = 0, with P = /C,f/. 

where a = /3/2, and /3 are defined just above. The solution of equation (16.121) is obtained in [2T] 
by reducing it to an obstacle problem for the s-Laplacian operator. On the other hand, Biler, 
Imbert and Karch [H] have obtained the explicit formula for a self-similar solution in the form 

(6.13) u{x, t- Ci) = t-"(Ci - h x'r^"/")^"^ 

with ki,Ci > and a = n/{n + 2 — 2s) as before. These solutions bear a resemblance to 
those of the standard porous medium equation [591 [ISl EO] • This is why they have been termed 
Barenblatt solutions for the FPME. 

Selfsimilar solution for s = 1: This approach works for s = 1 and gives a = 1, P = 1/n, 
where we replace ICg by /Ci = (— A)^^. We then recover the vortex patch solution (16. 8p . It is also 
easy to see that the family of self-similar solutions (I6.13P pass to the limit s — t- 1 to produce the 
vortex patch solution (16.81) (after careful selection of the constants). This is an explicit example 
of application of Theorem 13.11 

Comment. For the reader's convenience, we give more details in the case s = 1. The solution 
of the obstacle problem has to satisfy U = —AP = 2na = njH = 1 in the coincidence set \y\ < R, 
while AP = outside (for \y\ > R). Therefore, in the outer region 

p = D|y|2-"^ D>0. 

agreement of the two expressions aX y = r gives D = R^/n{n — 2) and R^ = 2{n — 2)C, and 
we recover the solution (16. Sp . 

6.4 Other examples of vortex patches 

We may next construct many solutions of the radial problem as rarefaction waves, which are not 
very different from the elementary vortex patch. We take an initial configuration formed by a 
ball patch at the origin plus an annular patch around, i.e., 

Uo{r) = Ci for < r < Pi, Uo{r) = C2 for R2 < r < P3, 
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with ci, C2 > and < i?i < i?2 < -R3, Wo('") = otherwise. Then 

M(r, 0) = (ci/n)r" for < r < i?i 

= {ci/n)BJl for Ri < r < R2 

= {c2/n)r^ + {cJn)Rl - {c2ln)Rl for R2 < r < R^ 

= {c2/n){Rl - i?^) + {c^/n)Rl for r > R:, 

Note that the continuity of M requires that we choose C2-R2 = cii?". At later times the solution 
keeps the same shape, i.e. the support of u{r,t) is still the union of a ball and an annulus (both 
expanding). u{r,t) and M{r,t) are explicit, in regions that are separated by interfaces: 
(i) The first region is an inner ball where u and M have the form 

Ci 1 Ci T""" r"' 

u,{r,t) = —^ = — ; Mi(r,t) 



Cit + 1 t + Ti' ^'^ ra(cit + l) r2(t + ri) 

(ri = 1/ci). The bounding circle (centered at the origin) has radius Si{t) satisfying S[(t) = 
viSiit),t) = M(^i(t),t)/^i(t)"-i with Si{0) = Ri, hence 

S,{t) = R^{cit + 1)1/" = i?iq/"(t + rl)l/^ 

(ii) There is then an annular region with M{r,t) = ciRi/n and u{r,t) = 0, bounded by this 
interface and the circle of radius 

S2it) = RMt + (i?2/i?l)")'/" = R2{c2t + 1)1/" = i?2C2/"(t + r2)l/", = (i?2/i?l )Vi . 

(iii) There is then an annular region of positive density (initially R2 < r < R^) given by 



C2t + 1 t + r2 ' n(c2t + 1) n{t + T2) 

where T2 = I/C2 > ti. The outer interface of this annulus is the circle of radius 5*3 (t) = 

i?3(c2t+l)l/" = i23C2/"(t + r2)l/". 

(iv) Finally, in the exterior region, complement of the ball of radius S-i{t), M[r,t) is constant M 
and u is zero. 

Conclusion. In this example there are two connected components for the density that expand 
out with time. The inner patch is unperturbed by the outer mass, while the outer annular patch 
is spreading out in a self-similar way. The density in both patches is constant in space but 
the density outside depends on the density inside, time, and the ratio R2/R1. Both densities 
approach 1/t in first approximation as t — )■ 00. Moreover, the gap between the two patches is 

d(t) = C[(t + r2)i/"-(t + ri)i/l 
which relatively goes to zero like 0(t~^"~i''/") as t — ?■ 00. 
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6.5 Lack of comparison principle 

This example allows to construct a counterexample to comparison of densities. Let us take for 
first solution ui{r,t) the above example with parameters < -Ri < r2 < -R3. For second solution 
U2{x,t) let us take the elementary vortex patch with center xq such that |xo| = (-R2 + -R3)/2, 
with height equal or less than C2 and initial radius R4 much less than {R3 — i?2)/2. In this way 
we have mi(x,0) > U2{x,0). However, the support of U2 expands from the initial ball Br^{xo) 
while the support of Ui has a "hole" that eventually will "cover" that ball at some time. Hence 
comparison of densities does not persist in time. 

Note that the lack of a comparison of densities principle was shown for the solutions of the 
FPME with s > 1/2 in [20] • The argument there was less explicit due to the lack of explicit 
solutions. 

Recall that for radial solutions on the other hand we have comparison of the mass functions 
M{r,t). A similar result for FPME is known in dimension n = 1 where the approach of inte- 
grating in space is viable, [15]. 



7 Solutions with compact support 

We now return to non necessarily radial solutions to (11.51) but impose the condition of compact 
support on the data. This short section is essentially contained in [12], where more details can 
be found. We have kept the proof since some of the arguments we explain here are used in what 
follows. 

It is known [2D] that the solutions of the FPME for < s < 1 have the property of finite 
propagation so that initial data with compact support produces solutions with the same property 
for every fixed t > 0. We can check that the calculation of the supersolution used to prove the 
property for s < 1 can be extended to s = 1. However, a comparison argument with explicit 
solutions is easier in this case, by using the selfsimilar solutions constructed in Section 16.21 It 
works as follows. 

Theorem 7.1. Given bounded initial data Uq with compact support in the ball of radius Rq, then 
for every t > the solution u{-,t) is supported in the ball of radius Ro{l + H'Uolloo^)"^'^"- 

Proof. We consider the characteristics that according to [1] transport the solution. Using the 
equation we see that the speed of such lines is given by 

(7.1) vixo,t) =Vp{xo,t) = {V{-Ar'u){xo,t) = C, [ ^^^L^^^^^-^l dy 

where Ci = l/{{n — 2)un), with a;„ the volume of the unit ball in M". So it is bounded at all x 
and t and we conclude that the initial support propagates only to a finite distance in any finite 
time interval. 

Take now t > fixed and assume that we are at Xq, a point on the free boundary or boundary 
of the support fl{t) = {x : u{x,t) > 0}. We consider an extreme point, i.e., such that |xo(t)| = 
sup{|x| : X G dQ(t)}. We may assume after rotation that xq = {R{t),0, ■ ■ ■ ,0). We calculate 
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the velocity along the radius direction (which is the velocity with which the boundary separates 
from the origin). Then, 

(7.2) vi{xo,t) = Ci / ■ dy, 

The integrand has a positive sign for every y G Q(t). Therefore, the worst case is when we 
replace u by its maximum estimate given by Theorem 15. H to get 

\v,{xo,t)\ < [ ^^rfy, 

+ ^ jBnio) Fo — y\ 

where r = 1/||mo||oo- Using scaling we get 

t + r Jb,{o) |ei -yr 

with ei = (1, 0, ■ ■ ■ ,0). Next, we notice that this is the integral that would be obtained for the 
self-similar solution of Proposition 16.11 After a calculation in radial coordinates we have 

which, integrated in time, yields the estimate R{t) < Ro (^)^^", i.e. the result. □ 

Remark 7.1. We have really proved the following: Let u he a hounded nonnegative weak solution. 
Assume the support of the initial data is contained in a hall of radius Rq. Let U he the vortex- 
patch solution with same initial L°° norm and initial support Br^-^^O). Then the support of u at 
time t is contained in the support of U at time t. 

We conclude the section with a similar statement for solutions which are not necessarily initially 
bounded. 

Proposition 7.1. For every initial measure fi G A^^(M") with compact support, there exists a 
weak energy solution with compact support at all times and satisfying a hound of the type (13. 3p . 

Proof. The existence of a solution satisfying (13. 3 p is proven in Theorem 13.11 To check it has 
compact support, let us approximate first the initial data by smooth initial data with same 
compact support and same bound M for their total mass J uq dx. Then, by conservation of 
mass, for every t > we have the two a priori estimates ||u(t)||j;^i < M, ||M(t)||oo < lA- By 
the Nash-Sobolev embeddings we then get Vp(-,t) G L°°(]R"). A simple scaling argument then 
shows that 

(7.3) ||Vp(-,t)|U < C\\Ap{-MT\\^v{:t)\t-'^''' < CiMVn-("-^)/" 

Using the description by characteristics (valid for the approximating solutions, as we saw in the 
previous theorem) and integrating, it follows that the supports of the approximating solutions 
(compact by the previous theorem) are included in a ball of radius 

(7.4) R{t) <Ro + CaM^/"^^/'^ , 

where C2 is independent of the approximation (and neither does Rq). The conclusion follows by 
passing to the limit. Note that we estimate the growth of the support like 0(t^/") as expected 
in the best case. □ 
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8 Asymptotic behaviour 

The class of nonnegative solutions with finite mass has a property of simple asymptotic behaviour 
of self-similar form that is typical of diffusion processes, and also of some conservation laws. We 
want to prove that the self-similar vortex patch solutions of Section 16.21 are attractors of the 
finite-mass solutions. This has been proven in [12] in the case of compactly supported initial 
data. We extend it here to solutions with finite second moment and to all radial solutions. 

8.1 Asymptotics by the entropy method 

By using the entropy method we can prove asymptotic behaviour for solutions with general data, 
with the only extra requirement that the second moment of the initial mass distribution has to 
be finite. The proof follows the outline of the one done for the diffusive process with s < 1 in 

m- 

In order to consider the question of large-time behaviour it is better to switch to the renor- 
malized fiow. More specifically, we introduce a change of the space, time and density variables 
corresponding to the expected large time behaviour, as follows: 

(8.1) y = x/{t + 1)1/", r = log(l + t), U{y, r) = (t + 1) u{x, t) , 
as well as the corresponding change for pressure and velocity: 

Piy, r) = (t + l)("-2)/>(a;, t), V{y, r) = (t + 1^'^-'^^^, t) = -VP{y, r). 

As a model for the asymptotic behaviour we consider the elementary vortex of Subsection 16.21 
with initial height 1, i.e. , 

u,ix,t) = -^ for xeBR{t), i?(r) = i?o (t + I)'/". 

and u^{x,t) zero otherwise. We further choose Rq so that the total mass M is the same: 

j u^{x,t) dx = J uo{x,t) dx, 

(i.e., oJnRo ~ ll^olli = M). The corresponding renormalized density is 

(8.2) U4y,T) = x\y\<Ro, 

and the corresponding pressure, defined by P^{r,y) = (— A)~if/*(r, y) satisfies 

P,(r,y) ={C- {yV2n)) and - VP.(r,y) = | in BrM- 

In the proofs below we will consider the solution U {y, r) of the renormalized fiow associated to 
a solution u{x,t) via (18. ip . Some direct calculations show that U solves 

(8.3) Ur = V -[UiVP + y/n)], AP + U = 0, 
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where V and A indicate now differentiation w.r.t. y. Note that this is a mass conservation law 

Ur + V ■ (Uv) = 0, v = - VP - y/n, 

where the velocity is the sum of the renormalized particle velocity and the confining velocity 
—y/n produced by the change of variables. 

We next state the asymptotic behaviour result. 

Theorem 8.1. Let n > 2 and assume that uq G L^{W'), uq > 0, and uq has a finite second 
moment, J^„y^uo{y)dy < oo. Then, the speed and density converge to equilibrium after nor- 
malization, as follows: 

(i) U stabilizes towards the elementary patch with the same mass, in the sense that 
(8.4) lim \\U{y,T)-U4y)\U= lim t) = . 

r— >oo t— >oo 

We also have strong convergence U{y,T) — ?■ f/*(y) in all the Lp(M"') spaces, 1 < p < oo, and 
weak-* convergence in L°°(M"). 

(a) The corresponding velocities converge 



15) lim V{y,T) ^ 



T— >00 



n 



= 0. 



Proof, (i) By displacing the origin of time by one unit, we may assume that our solution satisfies 
the a priori bound u{x, t) < hence r) < 1. The mass of the renormalized solution will 
be M at all times r > 0. This time displacement also allows to assume that the initial energy is 
finite, i.e., Vpo(l/) e L'^iW). 

The boundedness of the orbit r) in fl L°° implies that as r — )■ oo there is a sequence 
such that U{-,Tk) converges weakly to some Uoo{y), and also < Uoo{y) < 1. 

(ii) If n > 3, we consider the "entropy" 



(8.6) Ent{U):=l-[ {\V + -U y') dy 



which is the version for the renormalized flow of the energy (13. ip . Recall that VP = —V and 
V ■ V = U . Not also that under the stated assumptions j\WP\'^dy = JUPdy. It is easily 
calculated that the entropy evolves according to the rule 



^Ent(f/(r)) = -D{U{r)), D{U{r)) 





v-y- 




n 



2 

dy. 



(See [2T] for a similar entropy calculation for s < 1.) 

(ii') If n = 2, this entropy is not finite and its definition needs to be modified according to 
Section 13.21 we set Uq to be any smooth compactly supported nonnegative function such that 
jUo = JU, and let Pq = {-A)-^Uo. We then define 

Ent(f/) := £^ Q|V(P-Po)|2 + f/o(P-Po) + i^f/y2^ dy 
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which is finite. Computing as in Section 13.21 we find that 

^Ent(t/(r)) = -D{U{t)) 
still holds, and we can continue the proof in the same way as for n > 3. 

(iii) The fact that Ent(t/(r)) is nonincreasing along the orbit of the renormalized fiow implies 
that there exists the limit 

(8.7) E^U) = lim Ent(f/(r)) . 

T—^OO 

Moreover, the integral J U \ V — {y / n)\'^ dy dr is convergent as an integral in an infinite time 
interval. Another consequence is that the second moment JU{y,T)y'^dy is bounded for all r 
(by 2n Ent([/(0)), to be precise), and this control of the behaviour as y — )■ oo allows to make 
sure that in the limit / U^{y) dy = M. 

(iv) We now take limits as — )■ oo not only in the family {U{y,Tk)}k, k — oo, but in the 
family of orbits {U{y,T + Tk)}k, for |/ G M" and < r < T for fixed T > 0; A; — ?■ oo. The 
compactness (which can be simply reproduced from Theorem 13. ip allows to show that along a 
further subsequence that we also call ^ cx) we get 

U{y,T + n) -^U{y,T), 

in the weak sense, as well as V{y,T + r^) — )■ V{y,T) strongly. Notice that U{y,0) = Uooiy)- 
Since (18.71) holds, we also have 



Ent(f/(r)) = E^{U), a constant; and D{U{t)) = . 
The last expression means that for a.e r > and a.e. y we have 

U{V-iy/n)) = 0, 

which according to the renormalized evolution equation (18. 3p and V = —VP means that drU = 
0. Therefore, the limit orbit is a stationary solution of the renormalized equation. 

(v) We claim that our candidate U^{y) of (18. 2p is the only stationary solution satisfying the 
conditions J U{y) dy = M and < U < 1 and V{y) = y/n on the support of U. 

First of all the a priori bound for V in L°° (due to integral kernel estimates on the representation 
formula performed in Formula (17. 3p ) implies that V = y/n cannot be true for large y, therefore 
U must be compactly supported. Once we know the compact support property, it is already 
proven in [12] that the asymptotic limit is f/ = f/* (this also follows easily from our analysis of 
Theorem 17.11) . 



(vi) The convergence of U{-,t) towards f/* is at this point weak in L^. But we also have 
supyU{y,T) < 1, this unilateral bound allows to improve the weak convergence to strong con- 
vergence in the ball i?R„(0), i.e., \\U{y,T) — U^{y)\\i — )■ in that ball. But the mass of 
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contained in that ball is the whole mass, so that U {y, r) must converge to zero in the com- 
plement of the ball. This proves the main asymptotic formula (18. 4p . Since the solutions are 
uniformly bounded, the convergence in all the spaces also follows for 1 < p < oo. However, 
the convergence in is impossible due to the discontinuous form of the limit function. 

(vii) The convergence of the velocities now follows from properties of the convolution formula 
for V in terms of U. 

(viii) We still have to check that in displacing the origin of times from t = to t = 1 we have not 
lost the property of finite second moment. Indeed, using the equation and several integrations 
by parts, we have 



d 
di 



ux dx = / V ■ (nVp) X dx = —2 / x ■ [uWp] dx = 2 I x • VpAp dx 



-2jVp- V(x ■ Vp) = -2 / I Vj9|=^ -fx- V| Vp| 
= (n-2)/|Vp|^d. = (n-2)/.pdx. 
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If n = 2 then we are done. If n > 3, we have found that the time derivative of the second 
moment is a constant times the energy (13.11) which is itself nonincreasing in time, hence bounded 
independently of time. It follows that the second moment remains bounded for all times. □ 



8.2 Asymptotic behaviour for radial solutions 

We will perform next the proof that the large-time behaviour of finite-mass solutions is given by 
the self-similar vortex patch solutions of Section 1613) in the case of radially symmetric data, where 
the Burgers' reformulation allows to use familiar tools. We do not need any further restriction 
on the class of initial data. A reference for asymptotic convergence of conservation laws is |46] . 
Section 3 of that paper applies directly in the case n = 1 and it also does when n > 2 after 
replacing the space variable r by s. Since we also want to estimate u{r,t), which is a weighted 
derivative of the mass function M{r,t), we give the details of the argument. 

Theorem 8.2. For any radial solution to (11.51) with integrable initial data, letting M{x,t) be its 
"mass function" as in (16.11) and assuming ma,Xr>o M{r,0) = Mq, we haveO < M{rt^^"',t) < r"/n 
and 

(8.8) |r"/n-M(rti/'^,t)| ^0 

uniformly m G M and t — )■ oo. Moreover, if Uq is compactly supported, the set of x such that 
M{\x\,t) < Mq, i.e. the support ofu{x,t), is included in a disk centered at the origin of radius 
S{t) such that S{t) > (nMot)^/" and 



(8.9) 



S{t) - {nMoty/'' ^0 as t ^ oo. 
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Proof, (i) We first check the upper bound. First, since M{r,t) is increasing in r, it has a hmit 
huir-^oo = sup^ M(r, t) which is also equal to J u{r, t) r'^'^dr = Mq hence is constant in time. 

In view of the universal estimate of Section[5]we have u{x,t) < 1/t so that M{r,t) < r"'/nt and 
since limj.^.oo M(r, t) = Mq, we deduce in view of the above that S{t) > {nMotY^'^ (this includes 
that possibility that 5(^)00. This settles half of the estimates. 

(ii) In order to bound the solution from below we consider a vertically displaced function 

Mi{r,t) = M{r,t)+e, e > 0. 

Then, 

i.e. Ml is a supersolution of the original Burgers' equation. If Rq be determined by the condition 
M[Ro, 0) > Mq — e. Rq depends on e and may be large but is is finite. Then 

Mi(r, 0)>e for < r < Rq, Mi(r, 0) = Mq + e for r > i?o 

We now consider the vortex patch solution U2{r,t) with initial data ^XBr and total mass Mq. 
For that we must have R^/{nT) = Mq. We let M2{r,t) be its mass function. Then 

If we choose r > so that e > Rq/ut , then we will have Mi(r, 0) > M2(r, 0). By the maximum 
principle for (16. 2p . we deduce Mi{r,t) > M2(r, t) for all r > and all t > 0. This means that 

M(rtV«, t) > -e for r < [n(Mo + e){l + {r/t))]'/^ . 

n{l + (r/t)) 

Passing to the limit t — )■ 00, the term in r disappears, and letting then e we easily arrive at 
both estimates (18. 8p . 

(iii) The estimates on the support for compactly supported data are derived in Section 3 of [12] , 
even with rates of convergence. □ 

As a consequence of this result, we have estimates for the radial velocity v = M/r"'~^ that we 
resemble the velocity of the elementary vortex patch V{r,t) = r/nt. They are as follows: the 
velocity satisfies the estimates 

r 

v{r, t) < — for < r < R{t) 

Tvij 

and 

(8.10) - - t("-i)/"w(rti/", t) ^ 

n 

uniformly in r < S{t)/e/'' as t ^ 00. We can also give estimates for the density. 
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Corollary 8.1. The density satisfies the estimates 

(8.11) u{r,t)<j forO<r<R{t) 

and tu{rt^^'^,t) — )■ XBfl^(o)(^) in L°° (with weak-* convergence) as t oo, where Ri = (nMoY^'^. 

Proof. It is convenient to introduce the rescaled density u{r,T) = tu{rt^^'^,t). We know that u 
is bounded uniformly for large times by Theorem 15.11 and its support is uniformly bounded, by 
Theorem 18.21 It therefore converges in the weak-* topology of L°°{W^) to a function U with 
compact support. Due to the relation between u and M, the limit U must be related to the 
rescaled mass function of the elementary vortex patch by r"'~^U{r) = tM'{rt^^"') = M'{r). We 
conclude that U{r) is the characteristic function of the ball of radius Ri. Since the limit is 
unique the convergence does not depend on the subsequence — j- oo. □ 

Remark 8.1. The example of the radial solution with an expanding hall together with an expand- 
ing annulus shows that the convergence of u toward the equilibrium state U{r) = Xbm^/^{o){^) 
cannot happen in the uniform topology. 
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